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NOMENCLATURE 


The  following  symbols  ere  employed  in  the  derivations  of  this 
thesis. 

Symbol  Explanation 

A^  Total  number  of  atoms  of  a  particular  species  in  the  reaotanta. 

a^  The  number  of  atoms  of  a  particular  species  contained  in  a 

product  of  combustion. 

Sj  The  stoichiometric  coefficient  in  the  equilibrium  equation. 

(a/G)  Mole  fraction  of  species  a. 

(a/G)P  Partial  pressure  of  species  a. 

C  Centigrade  temperature  scale. 

CONV  Convergence. 

Cal  Calorie . 

Cp  The  value  of  specific  heat  measured  at  constant  pressure. 

Cv  The  value  of  specific  heat  measured  at  constant  volume. 

CH2  Cp  value  for  hydrogen  at  some  T. 

CH2A  Cp  value  for  hydrogen  at  T  ♦  100  K. 

F  Fahrenheit  temperature  scale. 

G, g  Gibb's  free  energy  function. 

G  Total  number  of  moles  when  used  in  the  equations  of  coabusticm. 

gmole  Gram  mole . 

H, h  Enthalpy. 

Hf  Heat  of  formation  (chemical). 

Hr  Heat  of  reaction  (ohemical). 

Enthalpy  of  hydrogen  at  T. 
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WOMWCUTORS  (continued) 


Sjubol  JjjsglsMtion 

H2A  Enthalpy  of  hydrogen  at  T  ♦  100  K« 

H2S  Enthalpy  of  hydrogen  at  298* 16  K. 

K  Kelvin  temperate!*  scale. 

Kcal  One  thousand  calorias. 

Kp  Equilibrium  constant  basad  on  partial  pressure, 

k  Isentropic  exponent  (Cp/Cv). 

L  bork  other  than  expansion  work. 

Log  Nsparian  logarithm. 

N  Molecular  weight, 

m  Mass. 

N  Nuaber  of  moles. 

n^  The  number  of  moles  of  a  product  of  combustion. 

O/P  Oxidixer/fuel  molar  ratio. 

P,p  Pressure. 

PS  Combustion  pressure  (atmosphere*). 

PK1  Logarithm  value  of  Kpl  at  soma  T. 

PK1A  Logarithm  value  of  Kpl  at  T  ♦  100  K. 

Q  Thermodynamic  heat. 

Ru  Universal  gas  constant  (1.98718  ^^al  ^). 

RS  Fractional  value  between  hundreds  of  decrees. 

S,s  Entropy. 

8  Molar  ratio  ( oxidise r/fual). 


NOMENCLATURE  (continued) 


Symbol  Explanation 

T, t  Temperature  degrees  absolute. 

Tc  Adiabatic  flame  temperature. 

TOLER  Tolerance  level  for  the  matrix  solution. 

U, u  Thermodynamic  enternal  energy. 

V, v  Volume. 

Error  function  of  substance  i. 

(3,7 )  (Row  3 >  Column  7)  of  the  matrix. 

(2tlli6)  (Reference  2j  Page  1L8) 


CHAPTER  I 


INTRODUCTION 

A  rocket  engine  is  8  device  for  converting  the  thermochemicel 
* 

energy  of  one  or  more  propellants  into  exhaust  jet  kinetic  energy. 

The  term  propellant  is  applied  to  any  material,  solid  or  liquid, 
consumed  in  the  rocket  motor  for  producing  high-pressure,  high- 
temperature  gases.  It  is  generally  assuzaed,  and  there  is  some  evi¬ 
dence  to  support  the  assumption,  that,  under  the  pressure  and  tenqpera- 
ture  occurring  in  a  rocket-motor  combustion  chamber  and  nozzle,  the 
chemical  reaction  takes  place  under  conditions  which  approach  chemical 
equilibrium.  Accordingly,  in  the  ensuing  discussions  the  thermo¬ 
dynamic  properties  of  the  gas  flowing  through  the  rocket  motor  are 
assumed  to  be  those  from  thermochemical  equilibrium  calculations. 
(1:1426)* 

The  scope  of  this  paper  will  be  limited  to  the  equilibrium 
calculations  of  combustion  products  and  the  resulting  adiabatic  flame 
temperature.  The  problem  of  attempting  to  determine  the  equilibrium 

constituents  as  the  gases  progress  through  the  rocket  nozzle  will  not 

'  ■ 

Jtn  attempted  afa  this  time  although  tills  authui  lealisee  the 'very 
pronouns ad  pmohlem-whiah  exists.  3  0  ^  ‘2. 

As  the  static  temperature  decreases  during  the  systematical 
progression  of  the  gases  through  the  nozzle,  the  molecules  of  the 
propellant  which  have  dissociated  due  to  the  high  temperature  in  the 
combustion  chamber  tend  to  reassociate.  In  a  basic  rocketzy  problem 
frozen  flow  is  assumed,  that  is,  the  equilibrium  combustion  of  the 


*  Numbers  in  parentheses  refer  to  items  in  the  bibliography 
and  the  related  page  number. 
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gases  in  the  combustion  chamber  is  unaltered  during  the  expansion 
process  in  the  nozzle.  (Iib36)  In  the  condition  of  equilibrium  flow 
the  equilibrium  conqposition  corresponding  to  the  local  conditions  of 
pressure  and  temperature  is  obtained  in  each  section  of  the  noisla. 

The  preceding  statements  correspond  to  the  two  extreme  cases  of  frosen 
or  equilibrium  flow.  In  the  case  of  real  flow,  at  a  temperature  often 
exceeding  3000  K  near  the  chamber,  reactions  of  recombination  take 
place.  (2tl79)  However,  the  residence  time  in  the  nozsle  is  very  small 
and  these  reactions  may  not  have  sufficient  time  for  completion. 

Equilibrium  flow,  therefore,  is  not  obtained  and  the  real  flow  is 
situated  at  some  intermediate  point  between  the  two  extreme  cases. 

Before  the  age  of  electronic  computers,  the  problem  of  calculating 
the  combustion  properties  of  a  rocket  engine  was  very  tedious  and 
almost  impossible  to  solve.  Since  the  present  day  engineers  have 
access  to  the  electronic  coiqputers,  it  is  only  fitting  that  a  thesis  be 
presented  in  a  form  illustrating  the  solution  of  combustion  problems 
by  use  of  the  electronic  computers. 

A/ 

tallQMiflg  nTHfirntr** am  As  a  solution  of  a  relatively  simple  .  ft 
combustion  problem  consisting  of  only  six  products  of  combustion  JBy 
-edapAy  understanding  and  following  the  illustrated  procedures  a 
problem  consisting  of  any  number  of  products  of  combustion  may  be 
computed  and  solved  to  any  desired  degree  of  accuracy. 

r 
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CHAPTER'  II 


HEAT  OF  REACTION  (AHr) 

The  energy  released  or  absorbed  In  a  chemical  reaction,  under  the 
condition  that  the  products  are  returned  to  the  initial  temperature 
of  the  reactants,  is  defined  as  the  Heat  of  Reaction.  If  the  reaction 
takes  place  at  constant  pressure,  the  heat  of  reaction  is  equal  to  the 
change  in  enthalpy  during  the  reaction,  i.e„,  (H)  products  -  (H) 
reactants.  Of  course,  these  enthalpy  values  must  include  the  energy 
associated  with  the  chemical  bonds  of  each  substance  as  well  as  its 
thermal  energy  which  is  due  to  the  random  motion  of  the  molecules,  (ljtl) 
The  heat  of  reaction  at  constant  pressure  is  related  to  the  heat 
of  reaction  at  constant  volume  by  the  following  thermodynamic 
relationships 

dq  ■  du  +  pdv  ■  dh  -  vdp  -  dHr  (1) 

The  internal  energy  (U)  for  an  ideal  gas  is  a  function  of  tempera¬ 
ture  only.  From  this  fact  it  is  possible  to  make  the  asauiqotion  that 
for  an  ideal  gas,  the  finite  difference  of  internal  energy  CAU)  at 
constant  pressure  is  equal  to  the  same  value  at  constant  volume,  which 
defines 


(OJ)p  *  «NJ)V 

(2) 

However 

(AHr)y  -AU  v**^° 

(3) 

and 

(AHr)p  -AH 

U) 

Since  by  definition 


H  -  0  ♦  PV 


(5) 


k 


than 

AH  -AO  ♦A(PV)  (6) 

on  rearranging 

AU  -Ah  -A(pv)  (7) 

Then  by  substitution 

<Wr)v  -  (Mr)p  -A(PV)  -  £U)p  -  ^U)y  (8) 

which  can  be  written  as 

(&Hr)  -  (&Hr)  -  (AN  )  RuT  (9) 

v  p  gas 

Standard  Heat  of  Formation  is  the  heat  of  reaction  when  a 
compound  is  formed  from  its  elements  in  their  standard  state.,  The 
standard  state  of  an  element  is  its  stable  form  at  a  temperature  of 
296.16  K  or  25  C  and  a  pressure  of  one  atmosphere.  The  standard  state 
is  sometimes  referred  to  as  being  0  C  and  one  atmosphere.  The  standard- 
state  temperature  is  really  artificial,  for,  if  the  reaction  were  to 
take  place  isotherraally,  it  would  require  a  perfect  heat  sink  which 
removes  heat  from  the  syatem  as  fast  as  it  is  liberated  by  the  chemical 
reaction. 

When  a  substance  is  combined  with  oxygen  the  heat  of  reaction  is 
called  the  heat  of  combustion  since  the  reaction  is  one  in  which  the 
oxygen  is  the  oxidizer.  The  formation  of  carbon  dioxide  from  its 
elements  in  their  standard  state  is  as  follows 

C  (graphite)  ♦  02(gas)  -  02(gas)  1  (AHf )C02(gas)  (10) 

Here  the  carbon  must  be  in  the  graphite  state  since  that  is  the 
stable  form  of  carbon  at  25  C,  likewise  the  oxygen  and  the  carbon 
dioxide  must  be  in  the  gaseous  state.  Under  these  conditions  the  heat 
of  reaction  is  the  standard  heat  of  formation  of  C02<  In  this  case 
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the  heat  of  reaction  can  also  ba  called  tha  heat  of  coabuation. 

The  sign  associated  with  tha  heat  of  raaetion  la  an  Indication  as 
to  whether  the  reaction  is  endothermic  (takes  up  heat)  or  exothermic 
(gives  off  heat).  In  this  country  normally  a  minus  or  negative  sign  is 
associated  with  an  exothermic  reaction.  This  convention  is  not  true 
all  over  the  world.  An  easy  way  to  determine  which  convention  is  used 
for  any  particular  set  of  tables  is  to  check  the  heat  of  formation  of 
water.  Water  will  always  have  an  exothermic  heat  of  formation  and 
from  this,  the  correct  sign  convention  can  be  determined  for  tha  remain¬ 
ing  substances. 


CHAPTER  III 


CHEMICAL  EQUILIBRIUM  AND  THE  EQUILIBRIUM  CONSTANT 

There  are  several  different  ways  of  defining  and  expressing  the 
equilibrium  constant.  For  example,  depending  upon  the  circumstances, 
the  equilibrium  constant  can  be  based  on  concentrations,  activities, 
fugacities,  or  partial  pressures  of  the  constituents  present  in  the 
equilibrium  mixture.  In  situations  where  the  perfect  gas  law  is 
closely  followed  by  the  constituents  of  the  mixture,  the  equilibrium 
constant  based  on  partial  pressures  is  the  most  useful,  i.e.,  the  most 
easily  applied. 

An  arbitrary  thermodynamical  system  which  is  separated  from  its 
surroundings  by  a  control  surface  which  is  a  perfect  thermal  insulator, 
will  now  be  considered.  It  is  desired  to  determine  from  the  thermo¬ 
dynamic  characteristics  of  the  system  whether  or  not  the  system  is  in 
equilibrium  or  undergoing  a  spontaneous  change  of  state.  As  a  first  step 
in  the  consideration  of  the  aforementioned  problem,  it  is  necessary  to 
determine  what  information  can  be  derived  by  applying  the  first  and  second 
laws  of  thermodynamics  to  the  system.  The  only  information  furnished  by 
the  first  law  is  that  one  type  of  energy  can  be  converted  into  another, 
and  the  ratio  of  the  conversion  is  constant.  The  first  law  gives  no 
information  regarding  the  direction  in  which  an  energy  transformation 
will  proceed.  As  far  as  the  first  law  is  concerned,  it  is  possible  for 
the  conversion  to  proceed  in  any  direction  that  is  permitted  by  the 
constraints  imposed  on  the  thermodynamic  system.  Thus  the  first  law 
gives  no  information  whatsoever  regarding  either  the  spontaneity  or  the 
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direction  of  a  physical  or  chemical  process.  ( 2  *529) 

The  second  lav  established  the  principle  that  it  is  impossible  to 
take  heat  from  a  system  and  convert  it  into  work  without  simultaneous 
changes  occurring  in  the  system  or  in  its  environment.  From  the  second 
law  it  is  learned  that  only  for  a  reversible  adiabatic  (isentropic) 
process  is  the  entropy  constant.  From  the  entropy  principle  it  is 
deduced  that  the  only  changes  of  state  (processes)  that  are  possible  are 
those  which  produce  an  increase  in  the  entropy  of  an  isolated  system.  The 
second  law  gives  no  information  indicating  whether  or  not  a  process  will 
take  place  if  the  constraints  imposed  on  the  system  do  allow  the  entropy 
of  the  system  to  increase.  Thus  it  is  seen  that  neither  the  first  nor  the 
second  law  provides  the  information  for  determining  (1)  whether  a 
thermodynamic  system  is  in  equilibrium,  and  (2)  if  a  spontaneous  change 
of  state  can  occur  in  a  thermodynamic  system. 

The  internal  energy  (U)  of  a  thermodynamic  system  is  analogous  to 
the  potential  energy  of  a  dynamic  system.  A  dynamical  system  which  is 
at  rest  can  be  set  into  motion  spontaneously  if  the  constraints  imposed 
on  it  permit  the  potential  energy  of  the  system  to  diminish.  The 
potential  energy  of  the  system  tends  to  convert  itself  epontaneously  into 
kinetic  energy.  Studies  show  that,  assuming  conditions  permit,  a 
dynamical  system  will  be  in  equilibrium  when  its  potential  energy  is  a 
minimum.  It  is  a  fundamental  principle  of  dynamics  that  a  dynamical 
system  tends  to  attain  the  state  where  ita  potential  energy  is  a 
minimum,  subject,  of  course,  to  the  geometrical  constraints  imposed  on 
the  system.  Thus,  the  internal  energy  (U)  of  a  thermodynamical  system, 
if  permitted  to  do  so,  will  be  trsnsformsd  into  other  forms  of  energy 
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having  a  nonpotential  energy  nature,  until  the  equilibrium  condition 
for  the  system  is  attained.  Hence,  a  spontaneous  change  of  state  can 
take  place  in  a  thermodynamic  system  only  if  it  occurs  in  ths  direction 
which  tends  to  make  the  internal  energy  of  the  system  a  minimum,  and 
the  system  will  not  be  in  equilibrium  until  the  minimum  value  of  U  is 
attained.  Hence,  one  fundamental  principle  that  gives  spontaneity  of 
a  process  is  that  the  process  must  proceed  in  the  direction  that  will 
result  in  the  U  of  the  system  proceeding  to  a  condition  of  U-min. 

From  the  second  law  the  change  of  state  for  an  isolated  thermo¬ 
dynamic  system  must  result  in  an  increase  in  the  entropy  of  the  system. 
Since  the  entropy  of  the  system  must  increase,  the  equilibrium  condition 
will  be  attained  only  when  no  further  increase  in  entropy  is  possible, 
that  is,  when  the  entropy  of  the  system  is  a  maximum,  subject  to  the 
constraints  imposed  on  it.  Hence,  a  natural  process  can  proceed  only 
in  the  direction  allowing  S  of  the  system  to  approach  S-max. 

From  the  preceding  it  followe  that  the  nece8sary  and  sufficient 
conditions  for  a  thermodynamical  system  to  be  in  equilibrium  aret  (1)  its 
internal  energy  must  be  a  minimum,  and  (2)  its  entropy  must  be  a 
maximum.  Hence,  a  spontaneous  process  la  possible  only  if  the  change 
of  state  can  proceed  ao  that  the  conditions  of  U->  U-min  and  S->S-max 
are  satisfied. 

The  free-energy  denoted  by  0  (Gibbs  Function)  is  defined  by 
G-H-TS-U-TS  +  PV  (11) 

Since  H,  U,  T,  and  S  are  thermodynamic  point  properties,  the  free- 
energy  function  G  is  likewise  a  point  property.  Hence,  a  finite  change, 
such  aeAG  depends  only  on  the  final  and  initial  states  of  the  qrstem. 
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The  exact  differential  dg  is  given  by 

dg  ■  du  -  tds  -  sdt  ♦  pdv  ♦  vdp  (12) 

Combining  the  first  and  the  second  lavs  of  thermodynamics  reaulta 
in  the  following  familiar  equations 

tds  •  dq  ■  du  ♦  pdv  ■  dh  -  vdp  (13) 

One  very  important  statement  which  is  quite  often  overlooked  by 
the  student  of  thermodynamics  is  that  these  equations  are  baaed  on  a 
pure  substance,  that  is,  a  substance  which  is  chemically  homogeneous  and 
fixed  in  chemical  compos ition.  If  the  state  of  the  system  is  such  that 
no  chemical  reaction  can  occur  without  an  external  stimulus,  then  the 
system  is  in  complete  (mechanical,  thermal,  and  chemical)  equilibrium,, 
hhen  in  such  a  state,  the  system  is  chemically  homogeneous  and 
invariant;  so  it  is  a  pure  substance  and  the  preceding  relationships 
will  apply.  However,  when  a  spontaneous  chemical  reaction  occurs  it  is 
an  irreversible  process  of  a  nonpure  substance  and  the  equations 
become  (6:U82) 

tds  >du  +  pdv  (H4) 

or 

tds  >dh  -  vdp  (15) 

In  the  statement  of  the  first  law  for  a  closed  system,  where 
du  -  dq  -  dw  (16) 

it  is  considered  that  dw  is  equal  to  pdv  *  dL,  and  defines  dL  as  work 
other  than  expansion  work;  the  tds  equations  for  a  chemical  reaction 
may  then  be  written  as 

tds  ■  du  ♦  pdv  *  dL  (17) 


or 
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tds  *  dh  -  Trip  ♦  dL  (16) 

The  free -energy  equation  is  then  eteted  as 

dg  •  du  ♦  pdv  ♦  rdp  -  tde  -  edt  (19) 

Then  substituting  for  tds  from  equation  (17) 

dg  *  du  ♦  pdv  ♦  vdp  -  (du  +  pdv  dL)  -  sdt  (20) 

which  results  in 

dg  *  vdp  -  sdt  -  dL  (21) 

In  the  special  case  where  dp  ■  0  and  dt  ■  0,  the  finite  change 
in  free-energy  for  such  a  process  is  denoted  by 

■  <°2  '  V»P  '  <22) 


This  equation  shows  that  the  change  in  the  free-energy  function 
(CkG)tp  represents  the  maximum  (nonexpansion)  useful  work  obtained 
from  a  constant  temperature  and  pressure  process.  It  also  shows 
that,  for  a  process  conducted  so  that  dp  ■  0  and  dt  •  0,  the  useful 
work^L  is  equal  to  the  decrease  in  the  value  of  the  free-energy 
function. 

The  free-energy  function  finds  its  principal  use  in  establishing 
the  criteria  for  predicting  the  equilibrium  of  chemical  and  physical 
processes.  For  a  system  that  is  in  equilibrium  at  a  given  temperature 
and  pressure,  the  value  of  the  free-energy  function  is  a  minimum. 
Moreover,  for  a  process  to  take  place  spontaneously  without  change  in 
temperature  and  pressure,  the  corresponding  change  in  the  free-energy 
function  must  be  negative.  Hence,  for  a  spontaneous  change  (with  dp  •  0 
and  dt  -  0) 


<°2  -  Vtp  <° 


(23) 


n 


The  foregoing  provides  a  criterion  for  the  possible  existence  of 
equilibrium.  If  ■  0,  a  spontaneous  change  of  state  is  impossible* 

The  foregoing  shove  that,  if  the  free-energy  function  at  a  second 
state  (state  2)  is  less  than  it  is  in  the  initial  state  (state  1),  the 
system  must  eventually  pass  to  state  2.  Thus  the  change  in  the  free- 
energy  function  is  a  measure  of  the  tendency  of  a  system  to  change  its 
state }  that  tendency  is  called  the  escaping  tendency. 

From  the  discussion  of  the  free-energy  function  the  following 
conclusions  may  be  drawn. 

I 

1.  The  only  chemical  reaction  that  can  proceed  spontaneously  are 
those  for  which  the  free-energy  function  decreases. 

2.  If  a  chemical  reaction  gives  no  change  in  the  free-energy 
function,  then  the  reaction  is  in  the  equilibrium  state. 

3.  Since  the  direction  of  a  chemical  reaction  depends  on  the 
concentrations  of  the  chemical  species  involved,  the  values  of  the  free- 
energy  function  for  the  different  species  are  related  to  their 

co ncentrations. 

It •  If  the  reactants  in  a  chemical  process  cannot  form  the  desired 
products  by  a  given  isothermal  process,  they  cannot  form  them  by  any 
other  process,  because  the  change  in  the  free-energy  function  is 
independent  of  the  path  since  it  ia  a  property. 

When  H2O  and  CO2  and  other  products  of  combustion  are  subject  to 
temperatures  of  1*300  F  and  higher  they  tend  to  disaociate  in  a  manner 
and  to  a  degree  determined  solely  by  the  temperature  of  the  mixture. 

When  the  producta  of  combustion  dissociate,  best  is  absorbed,  which 
decreaae  the  energy  available  to  accelerate  the  combustion  gases 
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through  the  rocket  engine.  Dieeocistion  is  a  reversible  reaction,  that 
is,  one  that  will  go  in  either  direction.  In  any  mixture  of  products  and 
reactants  both  dissociation  and  combustion  are  constantly  occurring* 

When  the  rate  of  dissociation  equals  the  rate  of  combustion  the  reaction 
is  in  chemical  equilibrium  and  an  equilibrium  constant  exists. 

In  the  general  chemical  reaction 

aA  ♦  bB^=i  cC  +  dD  (2h) 

A,  B,  C,  D  represent  various  substances  and  a,  b,  c,  d  represent  the 
corresponding  stoichiometric  coefficients.  ( 9 tLi98) 

For  this  example  it  will  be  assumed  that  all  the  substances 
involved  represent  ideal  gases.  The  double  arrow  indicates  that  when 
equilibrium  is  achieved  each  substance  will  be  present  in  a  certain 
amount,  and  each  substance  will  exist  at  its  partial  pressure.  As 
previously  stated  for  a  system  to  be  in  equilibrium  dL  must  equal  0, 
that  is,  no  work  is  being  done.  From  this,  equation  (21)  becomes 

dg  ■  vdp  -  sdt  (25) 

which  for  a  isothermal  process  reduces  to 

dg  •  vdp  -  dp  (26) 

P 

Then  upon  intergrating  from  p°.to  p  and  g°  to  g,  equation  (26) 
becomes 

g  -  g°  ■  RuT  log  R_  (27) 

P° 

By  letting  subscript  zero  refer  to  standard  conditions  at  one  atmosphere 

g  -  g°  •  RuT  log  p  (28) 

where  p  is  in  atmospheres. 

From  the  mixture  of  gases  as  given  in  the  previous  general 
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* 


equation  (2lt),  the  following  can  be  written 
gA  -  g°A  -  RuT  log  pA 

(29) 

gB  -  g°B  •  RuT  log  pB 

(30) 

g(,  -  g°c  -  RuT  log  pc 

(31) 

*  g°D  ’  *“*  108  PD 

(32) 

Nov  by  writing  an  expression  for  the  change  in  Oibb's  fret-energy 
function  as  the  reaction  goes  to  completion,  that  is,  starting  with 
reactants  A  and  B  and  ending  with  products  C  and  D  and  since 

AO  ■  G  products  -  0  reactants  (33) 

then 

AG  -  cgj,  ♦  dgQ  -  ag^  -  bg^  (3h) 

Substitution  in  the  pressure  equations  (29)  through  (32),  gives 
AG  ■  cg°  ♦  c  RuT  log  ♦  dg°  ♦  d  RuT  log 

-ag°  -  a  RuT  log  p  -  bg°  -  b  RuT  log  p  (35) 

*  a  B  B 

However,  AG°  is  defined  as 

AG°  -  cg°  ♦  dg°  -  ag°  -  bg°  (36) 

C  D  A  B 


and  equation  (35)  then  reduces  to 

c  d 


AG  -AG°  ♦  RuT  log  Pg. . 

Pa‘PB 


(37) 


When  the  reaction  is  in  equilibrium  (AG  •  0),  and  since AG°  is  a 
constant  and  RuT  is  a  constant,  the  quantity 


*-*r*#jB 


Ut 


log 


Pb 


(38) 


must  be  equal  to  a  constant  which  is  defined  as 

log  Kp  (39) 

This  results  in  the  following  equality 
-dG° 

Kp  -  e  m  (JlO) 

where  Kp  is  the  equilibrium  constant  expressed  by  the  partial  pressures 
of  the  combustion  gases. 

For  this  paper  no  attempt  will  be  made  to  consider  the  equilibrium 
constant  which  is  expressed  in  values  of  concentrations,  activities,  or 
fugacities. 


CHAPTER  IV 


CALCULATION  OF  THE  ADIABATIC  FLAME  TEMPERATURE 

All  methods  for  determination  of  the  adiabatic  flame  temperature 
and  the  composition  of  the  combustion  products  are  based  on  two 
principles:  the  First  Law  of  Thermodynamics  and  the  Law  of  Conservation 
of  Mass.  These  two  principles  when  used  in  conjunction  with  experi¬ 
mentally  determined  chemical  equilibrium  data,  yield  a  system  of 
simultaneous  equations  in  terms  of  the  gas  composition.  The  various 
methods  of  determining  the  gas  composition  differ  in  the  manner  in 
which  the  simultaneous  equations  are  set  up  and  solved. 

By  applying  the  First  Law  of  Thermodynamics  and  the  Law  of 
Conservation  of  Mass  to  a  chemical  reaction  it  is  possible  to  define 
the  term  Adiabatic  Flame  Temperature.  Assuming  that  the  chemical 
reaction  is  at  constant  pressure  and  adiabatic,  and  that  there  is  no 
shaft  work  performed,  the  First  Law  of  Thermodynamics  states  that  the 
enthalpy  of  the  reactants  at  their  temperature  prior  to  the  reaction 
(Tb)  is  equal  to  the  enthalpy  ofthe  products  of  the  reaction  at  the 
flame  temperature  (Tc).  Here  the  enthalpy  must  include  the  chemical 
energy  of  the  molecules  of  each  specie  as  well  as  their  energy  due  to 
random  motion,  since  changes  in  internal  energy  occur  as  a  result  of 
the  rearrangement  of  the  molecules  during  reaction. 

The  Law  of  Conservation  of  Mass  applied  to  a  chemical  reaction 
requires  that  the  number  of  atoms  of  each  element  remain  constant. 

For  exanple,  the  number  of  atoms  of  a  single  element  in  the  products  after 
the  reaction  must  be  equal  to  the  number  of  atoms  of  the  same  element 


16 


in  the  reactents  prior  to  the  reaction,  and  the  same  is  true  for  all 
other  elements.  Thus  in  accordance  with  the  definition  of  a  mole,  it  can 
be  said  that  the  number  of  atoms  of  any  particular  species  is  constant. 

With  the  aforementioned  laws  of  natural  response,  the  necessary 
equations  and  calculations  may  now  be  presented  for  solving  any  possible 
system  of  combustion.  For  this  combustion  problem  liquid  hydrogen  and 
liquid  oxygen  will  be  burned.  However,  by  simply  rearranging  some  of 
the  basic  steps,  aqy  oxidizer  -  fuel  combination,  may  be  used. 

Now  before  the  reaction  equations  can  be  written,  a  decision  must 
be  made  as  to  what  products  of  combustion  are  likely  to  be  present  in 
the  combustion  gas.  For  the  adiabatic  reaction,  the  equilibrium  flame 
temperature  will  be  very  high  and  a  considerable  amount  of  dissociation 
can  be  expected. 

A  basic  knowledge  of  chemistry  is  required  to  determine  the  possible 
products.  By  comparing  the  equilibrium  constants  of  the  selected 
products,  some  may  be  eliminated  by  observing  that  their  Kp  value  is 
several  orders  of  magnitude  less  than  the  values  listed  for  the  other 
reactions.  The  presence  of  other  combinations  may  be  discounted  when 
the  experimental  fact  is  considered  that  a  chemical  system  will  tend  to 
adjust  itself  so  that  the  constituents  are  in  their  lowest  energy  state. 
After  a  decision  has  been  reached  with  regard  to  the  products  of 
combustion,  the  reaction  equations  can  then  be  written. 

In  the  combustion  of  liquid  oxygen  and  liquid  hydrogen,  a  possible 
assumption  of  probable  products  are,  Hg,  °2»  H20>  0H»  an* 

h202.  03  and  H202 

can  be  eliminated  by  applying  the  above-mentioned 
rules  as  these  are  both  high  energy  forms. 
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The  following  equation 

1H2  ♦  s  02  ^  »  H?  ♦  b  02  ♦  e  +  d  0  ♦  e  H  ♦  f  OH  (ill) 

list  the  six  assumed  products  of  combustion.  The  value  of  s/l  will  give 
the  molar,  oxidizer  to  fuel  ratio.  This  can  easily  be  converted  to 
mass  ratio,  if  desired,  by  use  of  the  following  equation 

m  -  N  M  (U2) 

However,  for  this  problem  only  the  molar  ratio  will  be  considered. 

In  the  development  of  the  system  for  solving  for  the  molar  values 
of  the  products  of  combustion  it  is  necessary  to  regress  a  moment  to 
discuss  the  selected  system.  "The  errors  are  distributed  among  all  the 
species"  is  the  title  of  reference  ( 2 >H|2 )•  This  nethod  is  conqpletely 
general  and  can  be  used  irrespective  of  the  number  of  substances 
present  in  the  combustion  products.  The  system  is  well  suited  for 
solution  by  electronic  conpu-uers,  which  is  the  primary  reason  it  was 
selected. 

The  system  of  equations  used  to  determine  the  equilibrium  composition 
at  a  given  pressure  snd  temperature  is  formed  by  linear  equations 
representing  the  conservation  of  the  species  and  by  non-linear  equations 
which  can  be  written  in  logarithmic  form  and  represent  the  various 
equilibria. 

By  referring  to  equation  (lil)  and  utilizing  the  basic  equation  for 
the  balance  of  individual  species 

*1  “  2  \  d*3) 

the  first  three  equations  may  be  formed. 

(14*) 


1.  Equation  for  hydrogen  balance 
2«2a  +  2c  +  e  +  f 
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Sqmtlon  for  oxygen  balance 

2s«2b*c+d  +  f  (li5) 

3.  Equation  for  total  nolce 

Q»a  +  b*c*d*e*f  (l«6) 

For  any  combustion  problem,  the  number  of  products  of  combustion 
is  equal  to  I,  and  the  nunfcer  of  different  species  (atoms)  is  equal  to 
J.  Then  from  the  equation 

Y  «  I  -  J  U7) 

the  value  of  Y  will  determine  the  absolute  nunfeer  of  Independent 
equilibrium  equations,, 

For  this  problem 

Y  -  6  -  2  (i»8) 

or 

Y  -  k  049) 

Therefore,  the  number  of  independent  equilibrium  equations  is  h.  The 
basic  equation  for  equilibrium  in  logarithmic  form  is 

log  zi  ■  2  loS  n j  (50) 

ij  J  1 

which  allows  the  formation  of  the  following  four  additional  equations. 

The  values  of  Kp  listed  in  NACA  Report  1037  an)  also  in  the  Appendix 
are  based  on  the  constituents  being  formed  from  their  atoms.  Therefore, 
the  equations  become 

li.  HgO  Equilibrium  Equation,  which  is  formed  from 

2H  +  0  H20  (51) 


and  results  in 


Kp(H20) 


((e/0)P)«  (d/0)P 


(52) 
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which  by  rearranging  becomes 

Kp(H20)  -  fi* 
e2  d  Pz 

Then  by  equating 
Kp(H20)  -  Kpl 

In  logarithmic  form  equation  (53)  may  be  stated  as 

2  log  P  ♦  log  Kpl  ■  log  c  -  2  log  e  -  log  d  ♦  2  log  0 
5.  Up  Equilibrium  Equation,  which  is  formed  from 

2H^  H2 
and  results  in 


Kp(H2) 


(«/Q)P 

Tie/ ojp? 


which  by  rearranging  become  s 


(53) 

(510 

(55) 

(56) 

(57) 


kp(h2)  -  *  2  (58) 

vr  P 

Then  by  defining  Kp(Hg)  as  Kp2  and  writing  in  logarithmic  form, 
equation  (58)  may  be  stated  as 

log  P  ♦  log  Kp2  -  log  a  -  2  log  e  ♦  log  0  (59) 

6.  OH  Equilibrium  Equation,  Wiich  if  formed  from 
0  ♦  H  ^OH  (60) 

and  re suits  in 

Kp(OH) - -  (61) 

(d/0)P  (e/0)P 

which  by  rearranging  becomes 


Kp(OH) 


_f  G 
de  P 


(62) 
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Then  by  defining  Kp(OH)  as  Kp3  and  writing  in  logarithmic  form, 
equation  (61)  may  be  stated  as 

log  P  ♦  log  Kp3  ■  log  f  -  log  d  -  log  e  ♦  log  0  (63) 

7.  Oj  Equilibrium  Equation,  which  is  formed  fro* 

20  ^±r02  (61*) 

and  results  in 


Kp(02) 


(b/a)  P 

((d/G)P)2 


which  by  rearranging  becomes 


(65) 


Kp(o_)  -42  (66) 

*  d^  P 

Then  by  defining  Kp(Og)  as  Kpii  and  writing  in  logarithmic  form, 
equation  (65)  may  be  stated  as 

log  P  +  log  Kpli  -  log  b  -  2  log  d  +  log  G  (67) 

This  results  in  an  array  of  mathematical  logic  consisting  of  seven 
equations  and  seven  unknownso  Such  a  system  can  be  solved  by  trial  and 
error.  It  can  be  linearized  by  choosing  arbitrary  values  for  the  number 
of  moles  (n°)  and  then  solving  for  an  error  function  (X^).  The  error 
functions  are  used  in  the  following  equation 

njL  -  nj  (1  ♦  Xj)  (66) 


If  is  sufficiently  small,  equation  (68)  in  logarithmic  form  may  be 
written  a a 

log  ot  log'  n£  ♦  X°  (69) 

This  ia  true  by  considering  the  fact  that  the  log  (1  ♦  small  nuabar) 
is  approximately  equal  to  the  small  number. 

Equations  (1*3)*  (1*1*)*  (1*5)*  and  (1*6)  can  thus  be  written  in 
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the  fora 


.  _o  To 
'*i  ni  Xi 


(70) 


where 


A 


o 

1 


(71) 


Equations  (55) >  (5 9),  (63),  and  (67)  can  be  written  in  the  form 


I 


aJ 


log  z1  - 


log  n° 


where 

log 


log  n° 


(72) 


(73) 


The  transformation  of  the  variables  is  now  complete,.  The  new  system 
is  linear  with  respect  to  the  error  function  ( X ^)  and  can  thus  be  solved 
by  any  standard  method  of  solving  simultaneous  linear  equations.  The 
selected  method  is  well  suited  for  solution  by  electronic  computers. 

A  complete  description  for  the  selected  method  to  solve  simultane¬ 
ous  linear  equations  is  presented  in  the  Appendix.  The  program  is  a 
partial  replica  of  a  method  listed  in  the  IBM  listing  of  selected  programs 
for  the  1620  computer.  (5.0.007)  The  original  program  was  based  on  the 
original  FORTRAN  system  and  had  to  be  completely  rewritten  so  as  to  be 
suitable  for  use  with  the  present  FORTRAN  program. 

Table  I,  on  the  following  page,  consists  of  the  seven  linear 
equations  of  combustion  consolidated  into  a  matrix  fora.  A  detailed  , 
example  is  attached  to  Table  I.  It  is  available  for  additional  reference, 
as  required,  to  complete  the  formation  of  a  similar  system  of  equations 
necessai7  to  solve  other  combustion  problems. 
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The  value  of  X^  obtained  from  the  solution  of  the  matrix  provides 
the  new  value  of  which  is  computed,  for  the  second  approximation, 
from  the  equation 

ni1  -  ni°  (1  ♦  Xt)  (7b) 

A  very  important  consideration  is  that  the  first  approximation  may 
not  be  very  close  to  the  correct  value.  It  is  probable  that  duri(£  the 
first  few  trials  a  value  for  X^  of  less  than  -1.00  may  be  computed.  Hie 
program  will  not  accept  a  value  for  X^  of  -1.00  or  less;  a  value  greater 
than  -1.00  must  then  be  used.  In  this  problem  -0.995  was  utilized  for 
such  cases  with  complete  success.  After  the  assumed  values  are  within 
the  correct  range  a  value  of  less  than  -1.00  will  never  be  calculated  and 
the  calculated  values  of  X^  may  be  used  completely.  At  this  time  the 
computer  may  be  set  to  compute  continuously  until  the  desired  convergence 
is  obtained. 

When  the  calculated  error  functions  have  a  value  of  equal  to  or  less 
than  -1.00  and  are  then  used  in  equation  (7b),  a  negative  or  zero 
value  is  calculated  for  n^.  This  creates  errors  in  the  calculation 
resulting  from  the  computer  attempting  to  take  the  log  of  a  negative  or 
zero  number.  By  inserting  -0.995,  a  very  small  positive  value  of  n^ 
is  calculated.  This  answer  is  approximately  correct  since  it  would  be 
impossible  to  actually  have  a  negative  value  of  a  constituent. 

The  initial  selected  values  for  the  number  of  moles  of  the  products 
of  combustion  are  completely  arbitrary;  however,  if  any  information  on 
the  approximate  values  is  available,  it  should  be  utilized.  After  the 
problem  is  solved  one  time  it  is  advisable  to  use  the  calculated  values 


*5  tha  original  selected  values  for  the  number  of  mole  a  in  solving  the 
problem  for  the  new  temperature  selection  as  determined  by  Figure  U-l. 

This  will  allow  the  problem  to  converge  quicker  and  save  having  to 
replace  the  error  functions  that  are  equal  to  or  less  than  -1.00. 

Another  suggestion  is  to  require  the  sum  of  the  first  selected  values 
not  to  exceed  the  value  of  the  number  of  moles  of  the  reactants. 

The  convergence  of  the  system  is  given  by  the  expression 

CONV  -  /Ax  -  kx0/  ♦  / . /  ♦  /  In  !i/  ♦  / . 7  (75) 

Z1 

In  the  computer  program,  the  value  of  A^  -  A^°  is  set  equal  to  Al  and 
the  value  of  the  log  is  set  equal  to  Ah.  By  reference  to  the  initial 

equations  (Ui)  and  <s3.  it  is  quite  apparent  that  when  the  correct 
values  for  the  number  of  moles  of  the  products  have  been  confuted  Al,  A2, 
A3,  All,  A5,  A6,  and  A7  will  all  be  equal  to  zero. 

When  the  computer  does  not  have  the  ability  to  add  absolute  values, 
a  siirple  method  that  may  be  used  calls  for  adding  the  sum  of  ^  x*x  . 

The  system  is  solved  when  the  value  of  C  becomes  sufficiently  small. 

Hie  method  described  above  will  now  be  applied  in  the  solution  of  a 
problem  involving  the  selected  propellants  of  liquid  oxygen  and  liquid 
hydrogen. 

The  Kp  values  listed  in  NACA  1037  are  presented  in  logarithmic  form 
to  the  base  10.  The  logarithmic  program  on  the  1620  computer  is  based  on 
the  Naperian  Logarithm  Scale.  To  be  consistent  the  values  for  Kp  are 
multiplied  by  the  natural  logarithm  of  10. 

The  computer  program  is  subdivided  into  four  parts j 

1.  Solution  of  the  coefficients  and  the  values  for  the  solution  column 


of  the  matrix, 
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2.  Assignment  of  the  appropriate  values  to  columns  and 
rows  of  the  matrix  program. 

3.  Solution  of  the  matrix  program. 

U.  Solution  of  the  difference  in  enthalpies. 

The  value  for  the  total  number  of  moles  has  been  assigned  to  0 
and  GH  in  the  computer  program.  The  matrix  coefficient  (GH)  in 
position  (3,  7)  of  the  developed  matrix  created  some  confusion,  during 
the  early  stages  of  the  problem.  By  experimentation  it  was  found  that 
by  letting  GH  be  equal  to  the  new  value  for  the  total  number  of  moles 
calculated  from  the  value  of  XG  and  by  letting  the  G  equal  the  sum  of 
the  new  values  for  the  n^'s  the  problem  would  converge.  However,  by 
reversing  the  value  of  G  and  GH  the  problem  would  not  converge. 

The  computer  program  is  designed  to  automatically  replace  any 
undesired  error  function  with  -0.995.  This  ability  is  dependent  on  the 
position  of  Sense  Switch  3  viiich  sIbo  controls  the  type  out  of 
solutions,  When  Sense  Switch  3  is  in  the  off  position  all  values  are 
typed  out  by  the  typewriter.  During  the  type -out  of  the  solution  of  the 
error  functions,  if  one  has  a  value  of  -1.00  or  less,  it  is  changed  to 
-0.995.  A  message  that  this  step  has  been  completed  is  also  typed  out. 

When  the  operator  desires,  Sense  Switch  3  may  be  turned  on  to 
eliminate  typing  of  all  values  except  the  value  of  convergence.  This 
also  eliminates  the  ability  to  change  the  error  functions;  therefore,  it 
is  important  that  the  convergence  value  must  be  observed  to  be  decreasing 
before  this  step  is  performed.  If  Sense  Switch  3  is  turned  on  pre¬ 
maturely  and  an  undesired  value  of  an  error  function  is  computed,  an 
error  statement  will  be  typed  out  indicating  that  the  conqputer  has 


26 


b«en  requested  to  take  the  logarithm  of  a  negative  number.  The 
solution  must  be  restarted  ae  all  values  are  then  incorrect* 

As  the  convergence  approaches  the  desired  limit.  Sense  Switch  3 
may  be  turned  off.  This  allows  all  values  to  be  typed  out  which  are 
then  the  correct  solutions  for  the  number  of  moles  of  each  of  the 
products  of  combustion. 

By  turning  on  Sense  Switch  2,  the  computer  branches  to  the  fourth 
subdivision,  where  the  difference  in  enthalpies  is  computed. 

The  computer  also  has  been  programed  to  interpolate  linearly,  as 
necessary,  all  values  of  input  data.  The  value  of  RS  is  the  fractional 
value  between  hundreds  of  degrees.  If  it  is  desired  to  compute  for 
3528  K,  values  far  3500  K  and  3600  K  are  read  into  the  computer.  The 
value  of  RS  is  then  0.28. 

The  present  computer,  with  20,000  storage  locations,  did  not  have 
sufficient  locations  to  compute  the  value  of  the  molecular  weight  and 
the  isentropic  exponent.  A  second  program  was  written  to  perform  this 
computation.  The  last  step  of  the  first  program  is  the  command  to  punch 
out  on  data  cards,  the  solution  for  the  number  of  moles  of  each  of  the 
products.  These  solution  cards  are  then  used  in  the  second  program  with 
the  necessary  additional  data  to  compute  the  molecular  weight  and  the 
isentropic  exponent.  This  program  also  has  the  ability  to  interpolate 
linearly  the  necessary  data. 

For  the  design  calculations  in  this  problem  the  reactants  will  be 
assumed  to  be  at  standard  conditions  (298.16  K).  The  problem  is  than 
solved  for  a  particular  molar  ratio  and  a  selected  chamber  or  total 
pressure  (23  atmospheres).  The  necessary  values  are  taken  from 
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NACA  1037  by  assuming  a  particular  combustion  temperature .  Hie 
combustion  products  and  the  change  in  enthalpies  (AH)  are  then  calculated. 
AH  is  defined  as  the  enthalpy  of  the  products  taken  at  the  combustion 
temperature  minus  the  enthalpy  of  the  reactants  taken  at  the  temperature 
prior  to  combustion.  If  the  assumption  for  the  combustion  temperature  is 
correct,  AH  will  be  equal  to  zero.  If  AH  ^0  the  assumed  temperature  was 
too  high}  a  lower  value  must  be  chosen,  combustion  of  the  gas  must 
be  recalculated  for  the  new  temperature.  Of  course,  if  AH  <[0  then  a 
larger  value  for  the  temperature  must  be  selected.  Repetition  of  this 
procedure  two  or  three  times  will  usually  yield  the  correct  value  of 
temperature,  by  graphical  interpolation. 


Figure  li-l.  Interpolation  for  adiabatic  flame  temperature. 

If  the  first  two  assumed  values  of  temperature  give  values  of AH  on 
each  side  of  the  zero  line  as  shown  in  Figure  U-l,  an  approximate  value 
of  the  combustion  temperature  may  be  determined  by  linear  interpolation} 
however,  since  the  curve  of AH  vs  T  is  usually  concaved  up,  the  value  so 
determined  will  usually  be  less  than  the  actual  temperature. 

Linear  interpolation  for  the  proper  values  of  H  and  Kp  is  acceptable. 
It  should  be  noted  that  linear  interpolation  for  Kp  over  a  wide  tempera¬ 
ture  ranges  is  not  accurate,  and  logarithmic  interpolation  is  required 
for  adequate  approximations.  However,  the  equilibrium  constants  in 
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NACA  1037  are  tabulated  for  every  one  hundred  degrees  Kelvin  and  linear 
Interpolation  within  this  range  of  temperature  yields  sufficiently 
accurate  re suits . 

The  following  pages  consisting  of  Table  II  and  Table  III  have  the 
complete  conqputer  program  for  liquid  hydrogen  and  liquid  oxygen.  With 
the  previous  explanations  and  a  basic  knowledge  of  the  IBM  1620  FORTRAN 
system,  no  problem  should  exist  in  setting  up  the  program  to  compute  the 
combustion  properties  for  any  propellant  combination. 
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TABLE  II 

FIRST  COMPUTER  PROGRAM 


DIMENSION  2(7,8) 

81  PORMAT(//(F12»8)) 

32  FORMAT(//8HSOLUTION) 

33  FORMAT (LLHDIAGONAL  THAT  IS  LESS  THAN  OR  EQUAL  TO  TOLER) 

3h  FORMAT(2HIS  13) 

UO  FORMAT(ElU8) 

READiiO,  PK1,PK2,PK3,PK1*,PS,S, TOLER, RS 
READijO ,  PK1A ,  PK2A,  PK3A ,  PKL  A 
READl|0,A,B,C,D,E,F 
READ  LO,XA,XB,XC,XD,X£,XF,XG 
PK1-PK1  -  RS*(PK1  -  PK1A) 

PK2-PK2  -  RS*( PK2  -  PK2A) 

PK3-PK3  -  RS*(PK3  -  PK3A) 

PKL-PKL  -  RS*(PKli  -  PKliA) 

27  G-A+B+C+D+E+F 
A-A*(1„+XA) 

B-B*(1.+XB) 

C-C*(lo+XC) 

D*D*(l.+XD) 

E-E*(1.+XE) 

F-F*(1.+XF) 

GH-G*(1.+XG) 

G-A+B+C+B+E+F 

H-LDG(10.) 

Al-2.-(2.*A+2.*C+E+F) 

A2-2.*5-(2.*B+C+D+F) 

A3  -  G-( A+B+C+D+E+F ) 

AL-2.  *(  LOG( PS ) ) ♦H*PK1- ( LOG ( C ) -2 . *LOG ( E ) -L0G( D )+2 . *LOG ( G ) ) 
A5-(L0G(PS))*H*PK2-(L0G(A)-2.*LDG(E)*L0G(G)) 
A6-(LC>G(PS))*H*PK3-(LOG(F)-LOa(E)-LOG(D)+LOG(G)) 
A7"(LOG(PS))+H*PKii-(LOG(B)-2.*LOG(D)+LOG(G)) 

CON V-SQRT( A1*A1 )+SORT( A2*A2 ) +SQRT( A3*A3 ) +SQRT ( AL »Al ) +SQRT( A5*AS ) 
CONV-CONV+SQRT( A6*A6 )+SQRT( A7*A7 ) 

IF  (SENSE  SWITCH  3)50,51 
51  PRINT8l,Al,A2,A3,Ali,A5,A6,A7 
PRINT  81,A,B,C,D,E,F,G,GH 
50  PRINT  81,C0NV 

IF(SENSE  SWITCH  2)30,31 
31  X(l,l)-2„*A 
X(2,l)«0. 

X(3,l)-A 

X(l),l)-0. 

x(5,D-i. 

X(6,l)-0. 

X(7,l)-0. 

X(l,2)-0. 
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X(2,2)-2.*B 

X(3,2)-B 

XO»,2)-0. 

X(5,2)-0. 

X(6,2)-0. 

X(7,2)«l. 

X(1,3)-2.«C 

X(2,3)-C 

X(3,3)-C 

X(l4,3)«l. 

X(5,3)-0. 

X(6,3)-0. 

X(7,3)-0. 

X(l,li)-0. 

X(2,M»D 

X(3,U)-D 

X(l|,lj)—  1. 

x(5>lt)-o. 

X(6fli)— 1. 

X(7,l4)— 2o 

x(i,5)-e 

x(2,5)-0. 

X(3,5)-E 

x(Jj,5)-2. 

x(5,5)--2o 

X(6,5)-l. 

X(7,5)-Oo 

X(l,6)-F 

X(2,6)-F 

X(3,6)-F 

x(ij,6)-Oo 

x(5,6)-o. 

X(6,6)-lo 

X(7,6)-0. 

X(l,7)-0„ 

X(2,7)-Oo 

X(3»7)"-GH 

XU,7)-2. 

X(5,7)-l. 

X(6,7)-l. 

X(7,?)-l. 

X(1,8)-A1 

X(2,8)-A2 

X(3,8)-A3 

X(1i,8)-AJ4 

X(5,8)-A5 

X(6,8)«a6 

X(7,8)-A7 

N  -  7 

HI  »  N  ♦  1 
DO  U»  I  -  1,  N,  1 
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DIAQ  -  X(I,I) 

IF  (DIAQ)  3,  20,  k 
3  IF  (DIAQ  ♦  TOLER)  5,  19,  19 
14  IF  (DIAQ  -  TOI25R)  19,  19,  5 

5  DO  6  J  •  I,  N1 

6  X(I,J)  -  X(I,J)  /  DIAQ 
K  ■  1 

9  IF  (K-I)  11,  13,  11 

11  FELMT  -  X(K,I) 

DO  12  J  -  1,  Ml 

12  X(K,J)  -  X(K, J)  -  FELMT  *  X(I,J) 

13  x  ■  x  ♦  1 

IF  (K-N)  9,  9,  1U 
U4  CONTINUE 

IF  (SENSE  SViITCH  3)1,26 
28  J  -  N1 
TYPE  32 
DO  18  1-1,  N 
TYPE  I4O,  X(I,J) 

IF  (X(I,J)  ♦  1.0)  65,  65,  18 

65  X(I,J)  -  -.995 
TYPE  66 

66  FORMAT ( a 1H THE  ABOVE  VALUE  HAS  BEEN  CHANGED  TO  -.995) 

18  CONTINUE 
GO  TO  1 

19  TYPE  33 
TYPE  314,1 
GO  TO  5 

20  STOP 

1  XA-X(1,8) 

XB-X(2,8) 

XC-X(3,8) 

XD-X(1i,8) 

xe-x(5,8) 

XF-X(6,8) 

XG-X(7,8) 

GO  TO  27 

30  READ  L0,H2S,  02S,H2,02,H20,0,H,0H 
READ  l40,H2A,O2A,H20A,0A,HA,0HA 
OH  •  OH  ♦  RS*(OHA  -  OH) 

H  -  H  ♦  RS*(HA  -  H) 

0  -  0  +  RS*(OA  -  0) 

•»  H20  •  H20  ♦  RS*(H20A  -  H20) 

02  -  02  ♦  RS*(02A  -  02) 

H2  •  H2  ♦  RS*(H2A  -  H2) 

HTS-H2S+S«02S 

HTC-A«H2  ♦  B*02  ♦  C«H20  ♦  D*0  ♦  E*H  ♦  F«OH 
DELH-HTC-HTS 
PRINT  U1,HTC,HTS,DELH 
111  FORMAT  (//(Bill  •  8) ) 

PUNCH  a0,A,B,C,D,E,F 

STOP 

HID 
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TABLE  III 

SECOND  COMPUTER  PROGRAM 


67  FORMAT  (//(35HTHL  MOLECULAR  HEIGHT  OF  THE  MIX  IS  F10.5)) 

68  FORMAT (28HTHE  CP  VALUE  FOR  THE  MIX  IS  F10.5) 

69  FORMAT  (23HTOE  VALUE  FOR  CP/CV  IS  FIO.S) 

I4O  FORMAT(EU.8) 

READ  hO,A,B,C,D,E,F,S,RS 

READ  I4O,  CH2,  C02,  0120,  CO,  CH,  COH 

READ  hO,  CH2A,  C02A,  O120A,  COA,  CHA,  COHA 

G-A+B+C+D+E+F 

CH2  -  CH2  ♦  RS*(CH2A  -  CH2) 

C02  -  C02  ♦  RS*(C02A  -  C02) 

CH20  -  CH20  +  RS*(CH20A  -  CH20) 

CO  -  CO  ♦  RS*(COA  -  CO) 

CH  -  CH  ♦  RS*(CHA  -  CH) 

COH  -  COH  ♦  RS»(COHA  -  COH) 

CPM  -  CH2*A/G  >  C02*B/G  ♦  CH20*C/3  ♦  CO*D/G  +  CH*E/G  ♦  COH#F/G 
CVM  -  CPM  -  1.96718 
CK  -  CPM/CVM 

CMOL  -  2.*A/G  ♦  32.»E/G  ♦  18„«C/G  ♦  16.*D/G  ♦  1„*E/G  ♦  17.*F/0 

PRINT  67, CMOL 

PRINT  68, CPM 

PRINT  69, CK 

STOP 

END 
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CHAPTER  V 


RESULTS  OF  CALCULATIONS 

The  results  based  on  a  combustion  pressure  of  23  atmospheres  are 
the  primary  concern  of  this  report.  A  graph  relating  the  flame  temper¬ 
ature,  isentropic  exponent,  and  molecular  weight  of  the  products,  to  the 
oxidizer-fuel  ratio,  is  presented  in  Figure  5-2.  In  addition  a  very 
interesting  theory  of  rocketry  is  proven  -  the  temperature  of  the  gases 
depends  primarily  on  the  oxidizer-fuel  ratio,  and  to  a  small  extent  on 
combustion  pressure.  The  latter  depends,  of  course,  on  the  rate  at 
which  the  propellants  are  being  consumed  and  the  throat  area  of  the 
exhaust  nozzle.  The  curve  of  the  flame  temperatures  based  on  a  combus¬ 
tion  pressure  of  20  atmospheres  is  presented  for  demonstration  of  the 
effect  resulting  from  varying  the  combustion  pressure. 

The  program  is  set  up  for  the  IBM  1620  conputer  as  described  in 
Chapter  IV.  To  facilitate  the  calculations  of  the  initial  constants  and 
coefficients  for  the  matrix  program,  the  value  of  zero  must  be  read  in 
for  the  error  functions  XA,  XB,  XC,  XD,  XE,  XF,  and  XQ.  A  temperature  is 
selected  for  which  the  equilibrium  constants  are  taken  from  NACA  1037. 

In  addition  to  the  four  Kp  constants,  values  for  the  molar  ratio  (S), 
combustion  pressure  (PS),  tolerance  level  (TOLER),  and  a  value  for  the 
fractional  number  of  degrees  (RS)  have  to  be  read  in  to  the  computer. 

The  tolerance  level  and  its  use  are  explained  in  the  Appendix. 

Subdivision  three  of  the  first  conqouter  program,  consisting  of  the 
solution  of  the  seven  linear  equations,  has  almost  entirely  been  taken 
from  the  program  listed  in  the  Appendix.  The  input  section  has  been 
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deleted,  since  each  individual  position  of  the  matrix  has  been  identified 
in  subdivision  two.  The  portion  evaluating  a  zero  diagonal  has  been 
amended  because  the  only  way  to  adjust  the  rows  or  columns  of  the  matrix 
is  to  rewrite  the  program.  Care  should  be  taken  in  the  original  writing 
of  the  computer  program  to  prevent  a  zero  major  diagonal.  Specific 
instructions  are  given  in  the  Appendix. 

A  tolerance  level  value  of  0.00001  was  read  in  at  the  beginning  of 
the  program.  The  section  permitting  changing  of  the  tolerance  level 
was  omitted  from  the  computer  program.  The  tolerance  level  in  this 
program  is  not  very  important  since  the  calculated  values  of  the  error 
functions  are  only  used  to  create  new  assumed  selections  for  the  products 
of  combustion.  As  long  as  the  convergence  continues  to  decrease  the 
program  is  successful. 

The  computer  program  has  been  written  to  interpolate,  when  necessary, 
the  input  data  for  a  selected  temperature  falling  between  hundreds  of 
degrees.  An  example  may  be  illustrated  by  assuming  the  desired  tempera¬ 
ture  is  3528  K.  A  value  for  RS  of  0.28  must  be  entered  in  the 
computer  along  with  the  appropriate  input  data  selected  at  3600  K  and 
3500  K.  The  data  consisting  of  the  four  Kp  constants,  and  the  enthalpy 
of  the  constituents  are  taken  from  the  NACA  1037  tables  in  the  Appendix. 

The  values  listed  in  Table  IV  are  computed  from  an  oxidizer/fuel 
molar  ratio  of  0.5  and  with  a  combustion  pressure  of  23  atmospheres.  The 
computation  was  completed  for  temperatures  of  31*00  K,  3500  K,  and  3600  K. 
The  calculated  values  of  the  moles  of  the  products  are  given  in  Table  II 
along  with  the  value  of AH  computed  for  each  of  the  three  temperatures. 
These  values  of AH  are  plotted  on  the  graph  in  Figure  5-1,  and  then 
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graphically  interpolated  for  the  temperature  for  which  AH  is  zero.  The 
temperature  was  determined  to  be  3526  K.  The  problem  was  then  computed 
again  for  a  temperature  of  3526  K  and  the  resultingAH  was  within  the 
required  accuracy.  As  proof,  the  problem  was  then  confuted  at  3525  K 
to  get  an  estimate  of  the  change  of  AH  per  degree  Kelvin.  The  AH  at 
3526  K  was  +0.06076  while  the  Ah  at  3525  K  was  -0.10566.  This  was 
sufficient  proof  that  the  temperature  of  3526  K  was  within  the  experi¬ 
mental  accuracy  possible  in  the  combustion  chamber  of  a  rocket  engine. 


TABLE  IV 

RESULTS  OF  CALCULATIONS 
FOR  0/F  RATIO  OF  (0.5)  AND  23  ATMOSPHERES 


Temp.  K 

31400 

3500 

3525 

3528 

3600 

h2 

0.1325 

0.151)5 

0.1602 

0.1609 

0.1779 

°2 

.Oil  29 

.0ii91 

.0506 

.0508 

.055U 

h2o 

.6015 

.7629 

.7527 

.7515 

.7199 

0 

.0163 

.0257 

.0279 

.0281 

.03514 

H 

.0375 

.0518 

.0560 

.0565 

.070U 

OH 

.091)14 

.1132 

.1162 

.1188 

.1338 

G 

1.1271 

1.1572 

1.1656 

1.1666 

1.1928 

AH 

-6.5760 

-1.1)671 

-0.1057 

+0.0608 

+14.2998 

Since  the  second  computer  program  was  also  written  to  compute  a 
linear  interpolation  for  necessary  input  data,  it  is  only  required  to 
read  into  the  computer  Cp  values  for  the  six  products  of  combustion  at 
temperatures  of  3600  K  and  3500  K.  In  addition,  as  stated  before,  the 
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confuted  values  taken  from  the  first  program  for  the  number  of  moles, 
the  value  of  S,  and  the  value  of  RS  have  to  be  entered  as  input  data 
into  the  computer.  The  particular  values  of  molecular  might  were 
permanently  written  into  the  program. 

The  value  for  the  molecular  weight  of  the  mixture  is  computed  from 
the  following  equation 

M  -  (a/G)M  H2  ♦  (b/G)M  02  ♦  (c/G)M  HgO 

♦(d/G)M  0  ♦  (e/G)M  H  ♦  (f/G)M  OH  (76) 

This  value,  for  an  O/F  ratio  of  0.5  end  a  combustion  pressure  of  23 
atmospheres,  was  computed  as 

M  -  15-142976  (77) 

The  value  for  Cp  of  the  mixture  is  computed  from  the  following 
equation 

Cp  mixture  *  (a/0)Cp  H2  ♦  (b/G)Cp  02  ♦  (c/G)Cp  HgO 

♦(d/G)Cp  0  +  (e/G)Cp  H  ♦  (f/G)Cp  OH  (78) 

This  value  was  determined  to  be 

Cp  mixture  ■  ll.UjOSl  „  (79) 

gmole-K 


From  the  thermodynamic  identity  which  relates  the  specific  heat  at 
constant  valume  (Cv)  to  that  at  constant  pressure  (Cp),  Cv  will  now  be 
computed  from  the  equation 


Cv  ■  Cp  -  Ru 

(80) 

where 

*“  •  1-’ene  wn=K 

(81) 

Then  from  equation  (79)  and  (80) 

Cv  -  ll.lil*08l  -  1.98718 

(82) 
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or 

Cv  -  9.1)5363  — £2L_  (63) 

gmole-K 

The  isentropic  exponent  (k)  is  the  ratio  of  the  specific  heats 
and  is 

k  -  SE  (61)) 

Cv 


or 


k  .  11.1)1)061 

9.1)5363 


(85) 


then 


k  -  1.21020  (66) 

Table  V  consists  of  tabulated  values  from  the  calculations  at 
different  0/F  ratio  while  holding  the  combustion  pressure  constant  at 
23  atmospheres.  These  values  were  computed  with  the  same  procedures  as 
have  been  illustrated  for  an  0/F  ratio  of  0.5 

The  tabulated  results  of  the  change  in  adiabatic  flame  temperature, 
caused  by  changing  the  combustion  pressure  to  20  atmospheres,  is  present¬ 
ed  in  Table  VI.  Only  the  variation  of  the  flame  temperature  is  avail¬ 
able  in  Table  VI  as  this  is  the  only  variation  of  importance. 
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TABLE  7 

FINAL  RESULTS  OF  CALCULATIONS 
COMBUSTION  PRESSURE  23  ATMOSPHERES 


0/F  Molar  Ratio 

0.2 

0.3 

fmm 

0.6 

0.7 

*  Tc 

:  K 

2777 

3300 

3*)92 

3528 

31)95 

31)1)0 

h2 

0.5950 

0.391)2 

0.21)76 

0.1609 

0.1099 

0.0788 

-fJ 

o 

3 

°2 

.0000 

.0013 

.0139 

.0508 

.1106 

.1865 

O 

u 

a. 

h2o 

.3979 

.5679 

.6830 

.7515 

.7961 

.8270 

c 

•rt 

(0 

0 

.0001 

.0023 

.0131 

.0281 

.0395 

.0660 

a> 

rH 

O 

H 

.0121 

.01)86 

.0627 

.0565 

.01)1)8 

.031)5 

OH 

.0020 

.0272 

.0761 

.1188 

.11)32 

.151)1 

G 

1.0071 

1.01)15 

1.0961) 

1.1666 

1.21)1)1 

1.3269 

M 

8.31)05 

11.1377 

13.1)98  5 

15.1)298 

17.0399 

18.3909 

ad  Op 

10.2839 

11.0239 

11.3381) 

11.1)1)08 

11.1)600 

11.1)360 

rH 

o> 

CO 

° 

'o  Cv 

B 

8.2967 

9.0367 

9.3512 

9.1)536 

9.1)728 

9.1)1)88 

to 

k 

1.2395 

1.2199 

1.2125 

1.2102 

1.2098 

1.2103 

TABLE  71 

FINAL  RESULTS  OF  CALCULATIONS 
COMBUSTION  PRESSURE  20  ATMOSPHERES 


0/F  Molar  Ratio  0.2  0.3  O.U  0.5  0.6  0.7 


Tc  K 


2758  3280  31)70  3505  31)76  31)11 


■mmbmhHI 
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Figura  5-1*  Oraphical  Interpolation  for  Coaboatioa 
Flan  Taaparatova 


Ficon  5-t.  OrqphlMl  Basolta  For  Liquid  Hydrofoo 
•ai  liquid  Oxjrcoe 


JOLBCUUR  MSIGHT  (X) 


la 


TABLE  VII 

RESULTS  FROM  FIRST  PROGRAM 


-.30000000  A1 
-.lOOOOOOO  A2 
.00000000  A3 
l.lj 81*35890  All 
2.05136710  A5 
1.31*662590  A6 
1.14823U660  A7 

.1*0000000  h2 
.10000000  o2 
.1*0000000  h2o 
.10000000  0 
.30000000  H 
.1*0000000  OH 
1.70000000  G 
1.70000000  GH 


The  following  four  pages  art  the 
computer  results  from  the  following 
input  data* 


Initial  assumption 
for  products  of 
combustion. 


T  3528  K 
0/F  0.5 

TOLER  0.00001 

PS  23  Atmospheres 


9.761*69830 


Initial  convergence  value 


SOLUTION 

—523713UOE-00  XA 
- .  It  36OOOOOE-OO  XB 
.  11601*  7l*l*E*01  XC 
-.ll6951<lt5E*01  XD 
THE  ABOVE  VALUE  HAS 
-.U*979106e+01  XE 
THE  ABOVE  VALUE  HAS 
-.90008880E-00  XF 
-.1*2071*01*2E-00  XO 


Results  of  initial  calculations 
for  the  error  functions. 

BEEN  CHANGED  TO  -.995 
BEEN  CHANGED  TO  -.995 


-.15087310 

-.0171*51i  20 
.00000000 
•11.1*01*51800 
-7.1*1533000 
-6.55833300 

-8.15338600 


.190511*61* 

.0561*0000 

.861*18976 

.00050000 

.00150000 

.039961*1*8 

1.15306880 

.9SU7U132 


33.69989100 

SOLUTION 
-.18793932E-00 
.  18014  2300E-02 
-.75171W40E-01 
.10532358E+01 
.35893352E*01 
. 11329576E+01 
-.1*8720557E-01 

-.00000020 

-.00000000 

.00000000 

-6.573871*00 

-14.11718150 

-14.12958100 

-lj.  87263250 

.1514709145 

.05650175 

.7992271*1 

.00252661 

.006881400 

.085214251* 

1.10509170 

1.09689070 

19.69326700 

SOLUTION 
.314152510E-01 
-.73875990E-01 
-.314020070E-01 
.2i*102371E+01 
•20865257E+01 
.3I45I46I4I4OE-OO 
. 21717 281E-01 

-.00000010 

-.00000010 

.00000000 

-3.10101*200 

-1.91799860 

-2.0938191*0 

-2.36361*870 

.15999316 

.05232763 

.77203769 

.00861636 

.02121*765 
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.  111*69080 
1.12691320 
1.12909120 

9.1*7650880 

3.06712780  Here  Sense  Switch  3  was  turned  on. 

.1*1*1*07657 

.0111*7719 

.00000870  Here  Sense  Switch  3  was  turned  off. 

SOCtJTION 

.385711j00E-06  The  values  of  the  error 

-.58173300E-06 

.l*36l9620E-07  functions  are  now  very  small 

.15603613E-05 

.l*9l*081*66E-06  numbers. 

-.U*800950E-06 

.102l*5536E-06 

.00000010  Al 
-.00000000  A2 
.00000000  A3 
-.00000020  Al* 

-.00000010  a5 

-.00000020  A6 
-.00000020  A7 

.16086923  Ho  Final  values  for  the 

.05080107  02 

.7511*9282  HoO  products  of  combustion. 

.02813721  0 
.05650815  H 
.11876787  OH 
1.16657620  0 
1.16657620  (Hi 

.00000080  Final  convergence  value 

. 71556908E+02  Enthalpy  of  products  at  3528  K 

.711*96150B*02  Enthalpy  of  reactants  at  298.16  K 

.60758000E-01  AH  . 

STOP 


TABLE  VIII 


RESULTS  FROM  SECOND  PROGRAM 


THE  MOLECULAR  WEIGHT  OF  THE  MIX  IS  15. L 2976 

THE  CP  VALUE  FOR  THE  MIX  IS  11.M081 
THE  VALUE  FOR  CP/CV  IS  1,21020 
STOP 
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SOLUTION  OF  SIMULTANEOUS  LINEAR  EQUATIONS 


Introduction 


This  program  is  a  replica  of  program  number  (5*0.007)  from  the 
IBM  1620  Oeneral  Program  Library.  It  has  been  rewritten  as  necessary 
to  be  compatible  with  the  present  1620  (20K)  computer  and  the  FORTRAN 
programing  system  with  format  statements.  Although  this  program  was 
tested  by  its  author  using  the  sample  problems,  no  warranty,  expressed 
or  implied,  is  made  by  the  author  or  IBM. 

Description  of  the  Program 

A  knowledge  of  the  IBM  1620  FORTRAN  system  is  required  to  under¬ 
stand  the  following  program. 

This  program  solves  sets  of  non-homogeneous  simultaneous  linear 
equations  and  provides  either  typewriter  or  punched  card  output.  Up  to 
25  equations  in  25  unknowns  may  be  solved.  The  first  card  read  for 
each  set  of  equations  states  the  number  of  equations  following.  The 
program  automatically  sizes  itself  for  the  proper  matrix  dimensions. 

Any  number  of  sets  of  equations  may  be  solved  one  after  the  other  without 
reloading  the  program  deck,  khen  a  solution  has  been  obtained  the  word 
SOLUTION  is  typed.  Under  this  heading  the  solution  variables  are  typed 
one  per  line  in  the  order  in  which  they  appear  from  left  to  right  in  the 
original  equations.  Punched  card  output  is  identical,  with  one  card 
punched  for  each  line  typed. 

If  a  set  of  equations  is  submitted  which  has  more  than  one  solution, 
a  valid  solution  will  be  typed  or  punched,  but  no  indication  will  be 
given  that  there  are  other  solutions  as  well. 

If  the  equations  in  the  set  are  not  all  independent,  a  typewriter 
message  will  indicate  that  a  diagonal  element  is  zero.  No  solution  is 
developed  in  this  case,  as  there  are  infinite  solutions  to  the  set  of 
equations. 

Additional  information  regarding  the  characteristics  of  the 
program  is  contained  in  the  discussion  of  the  method  of  confutation. 

Method  of  Computation 

Oeneral  Description.  The  Jordan  method  of  elimination^  is  used  to 
accomplish  the  solution!  Briefly  this  method  consists  of  diagonalizing 
the  augmented  coefficient  matrix.  After  diagonalisation  the  right  hand 
or  augmented  column  contains  the  solution  values  for  the  variables. 

The  reduction  of  the  matrix  tikes  place  in  the  same  storage 
locations  where  the  original  coefficients  are  stored.  Thus  at  the 
completion  of  the  solution  the  original  coefficients  are  lost,  and  in 

1 Alston  S.  Householder,  Principles  of  Numerical  Analysis  (New  lorki 
McGraw-Hill  Book  Co., Inc.,  1953 ) >  PP*  68-72. 
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their  piece  remain  an  identity  matrix  and  the  column  of  answers. 

Division  by  a  Zero  Diagonal.  In  the  process  of  disgonalizstion 
the  elements  or  each  row  to  the” right  of  the  diagonal  and  the  diagonal 
itself  are  divided  by  the  diagonal  element.  The  diagonal  elements  are 
used  in  the  same  order  in  which  they  appear  in  the  original  matrix. 

No  attempt  is  made  to  select  these  pivotal  elements  by  magnitude.  Hence 
a  zero  must  not  appear  on  the  diagonal  of  the  matrix.  If  such  a 
condition  exists,  a  message  vdll  be  typed  indicating  which  diagonal 
element  is  zero  and  control  passes  to  reading  the  next  problem. 

The  possibility  of  zero  diagonal  elements  may  be  minimized  by 
observing  two  rules  in  laying  out  the  original  matrix. 

1.  The  upper  left  hand  element  should  not  be  zero. 

2.  No  diagonal  element  may  be  zero  if  all  elements 
in  the  row  to  the  left  of  it  are  zero. 

There  is  a  slight  chance  of  generating  a  zero  diagonal  element  in 
the  course  of  calculations.  Equations  which  are  not  independent  (a 
unique  solution  undefined)  will  also  generate  a  zero  diagonal.  If  it  is 
known  that  the  coefficient  matrix  is  non-singular,  that  is,  that  there 
is  a  unique  solution,  division  by  zero  may  be  easily  corrected  and  the 
problem  submitted  again.  A  simple  rearrangement  of  the  columns  will 
usually  suffice.  In  rearranging  the  columns  the  equality  vector  always 
must  be  kept  on  the  right.  As  indicated  later,  rearrangement  of 
columns  may  be  acconpliahed  simply  by  rearranging  input  cards.  No 
additional  punching  is  required. 

Tolerance  Warning  Level,  bhen  the  diagonal  element  is  very  small 
relative  to  the  values  to  its  right  in  the  row,  a  loss  of  accuracy  will 
result  from  round-off  errore.  In  this  program  the  operator  has  control 
over  this  situation.  Vlhen  the  absolute  value  of  the  divisor  is  less 
than  or  equal  to  some  tolerance  level  selected  by  the  operator,  a  warn¬ 
ing  message  ia  typed.  The  operator  then  has  the  choice  of  continuing 
with  the  solution  or  proceeding  to  the  next  problem.  If  he  should  elect 
to  continue,  he  does  so  with  the  knowledge  that  the  answer  may  contain 
significant  rounding  errors.  Of  course,  checking  the  answers  will 
quickly  show  the  accuracy  of  the  solution. 

Input  Card  Format 

The  first  card  of  the  input  for  each  set  of  equations  must  be  a 
card  with  the  number  of  equations  punched  according  to  the  FORMAT 
STATEMENT  (12). 

The  matrix  of  coefficients  ia  punched  column-wiee  beginning  with 
the  left-most  column  and  ending  with  the  equality  column.  Each  column 
must  begin  with  a  new  card.  A  card  must  contain  five  values  punched 
successively  according  to  the  FORMAT  STATEMENT  (Ell*. 8).  It  is  obvious 

•V 


w  rtjtunwy* 
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that  the  five  values  will  take  up  the  first  70  spaces  on  the  card*  The 
column  say  be  continued  over  as  many  cards  as  necessary.  (Five  cards  la 
the  maxinun  number,  as  the  program  is  limited  to  25  unknowns.)  If  the 
number  of  rows  is  not  a  multiple  of  five,  dummy  aero  values  muat  be 
punched  to  finish  out  five  entries  on  the  last  card  of  the  column.  Zt 
is  never  necessary  to  punch  more  than  four  dummy  sero  entries  for  a 
given  column.  The  number  of  cards  per  column  will  be  the  same  for  all 
columns  within  a  given  problem,  but  may  vary  from  problem  to  problem. 

The  dummy  entries  used  to  fill  out  the  five  entries  on  the  last  card  of 
a  column  need  not  be  zero.  These  entries  are  accessed  during  the  read 
in,  but  are  never  .accessed  during  calculation.  Hence  they  muat  be 
present,  but  their  value  is  immaterial. 


Output  Card  and  Typewriter  Format 


The  first  statement  for  each  problem  will  be  the  word  SOLUTION 
followed  by  the  answers  typed  according  to  FORMAT  STATEMENT  (Slit .8) .  If 
it  is  desired  to  have  the  answers  punched  with  the  same  format  state¬ 
ment,  Sense  Switch  1  should  be  turned  on.  The  word  SOLUTION  will  not 
be  punched;  however,  the  answers  including  the  word  SOLUTION  will 
continue  to  be  typed.  The  cards  will  contain  the  solution,  one  per  card 
in  the  order  in  which  the  variables  appeared  in  the  original  matrix. 

Sense  Switch  Settings 

Sense  Switch 


1 

2 

3 


ON  Solution  is  typed  and  punched. 

OFF  Solution  is  typed. 

ON  Next  problem  after  tolerance  stop. 

OFF  Continue  after  tolerance  stop. 

ON  Tolerance  level  entered  for  each  problem. 
OFF  Tolerance  level  entered  only  once. 
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ON  To  correct  error  in  typing  tolerance  level. 
OFF  Tolerance  level  entered  correctly. 


Normal  Loading  Procedure 


1.  Clear  Storage 

2.  Depress  RESET 

3.  Depress  LOAD  button 

U.  Depress  READER  START  when  hopper  is  empty 

5.  LOAD  DATA  is  typed  by  typewriter 

6.  Enter  data  cards 

7.  Depress  COMPUTER  START 

8.  Depress  READER  START 


Special  Loading  Instructions 


At  any  time  the  program  may  be  initialized  and  started  by 
branching  to  08300. 

Other  Instructions  and  Remarks 


Tolerance  Level,  khen  the  program  begins,  the  typewriter  will 
request  the  operator  to  enter  the  tolerance  level  as  a  floating  point 
number  in  the  form  (+.000000). 

If  the  entry  is  typed  correctly,  set  Sense  Switch  l\  OFF,  press 
RELEASE  and  COMPUTER  START.  The  data  cards  will  then  be  read. 

If  the  entry  is  typed  incorrectly,  turn  Sense  Switch  li  ON,  press 
RELEASE  and  START.  The  program  will  immediately  return  control  to 
the  typewriter  so  the  entry  may  be  made  again.  If  the  retry  is 
correct,  follow  the  procedure  for  a  correct  entry. 

If  it  is  desired  to  enter  a  tolerance  level  for  each  set  of 
equations,  turn  Sense  Switch  3  ON.  If  not,  turn  Switch  3  OFF,  and  the 
last  tolerance  level  entered  will  automatically  be  applied  to  all 
succeeding  sets  of  equations. 

khen  a  tolerance  warning  stop  occurs,  directions  for  the  use  of 
Switch  2  are  typed.  Only  at  this  time  is  the  setting  of  Switch  2 
interrogated. 

All  programed  stops  are  accompanied  by  typewriter  messages 
which  are  self  explanatory  and  which  indicate  the  required  sction. 
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SAMPLE  PROBLEMS 

The  following  four  problems  were  used  one  after  the  other. 

8*!  -5x2  *?x3  -  29 

+2x2  -2x3  -  -2 

-  xi  +9x2  -6x3  -  1 

Answers  x^  -2,  x2  -  3*  x3  -  Is 

X1  +2x2  +3x3  +5x^  -6x6  +7x?  -8xg  -  -28 

7x2  -2xj  ♦  x^  ♦  Xg  »  19 

2xx  +3x2  -Lx3  -5x^  +6x6  -  12 

X1  *  x2  4  x3  ♦  xu  +  x5  +  x6  +  x?  ^  xB  -  36 

10x^  +3x2  4jx3  +  xg  +2xg  -9x„  •  -is 2 

3x1  -3x2  -2x^  +2x^  +  xfi  -5 

-8x1  -9x3  +7x0  +6xy  -3xg  -  18 

+3x^  -5xg  -  -19 
Answers:  x^  ■  1,  x2  -  2,  x3  ■  3»  x^  ■  U 

x5  -  5,  x6  -  6,  x?  -  7,  xg  -  8 

3x^  +8x2  +6x3  +10XJJ  +L  2x0  •  -20 
+2x2  +  x^  +  5x0  -  i»0 

♦  lix^  ♦  6x^  -  -18 
2xx  +lsx2  +7x3  ♦  9X0  -  12 

X1  +^x2  *6x3  ‘  2xl  ■  7 

This  problem  will  cause  diagonal  3  to  be  zero.  It  was  included 
merely  to  illustrate  the  typewriter  message. 
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3xi  +8x2  +6  Xj  ♦lOx^  *U2x^  -  -20 

+2*2  *  +  5x5  -  1*0 

+  .Olxj  ♦  l*x^  ♦  6x^  ■  -18 

2xx  +1*X2  +7  x3  ♦  9xj  -  12 

X1  +Ux2  +6  x3  -  2x^  +10x^  -  7 

This  problem  will  make  diagonal  3  go  below  a  tolerance  level  of 
0.1.  It  was  included  merely  to  illustrate  this  typewriter  message. 


COMMENTS  ON  THE  TYPEWRITER  LOO 
FOR  SAMPLE  PROBLEMS 


On  the  typewriter  log,  which  follows,  several  points  should  be 
noted.  Sense  Switch  3  was  on  at  the  beginning.  Note  that  a  toler¬ 
ance  level  was  entered  for  each  of  the  first  three  problems.  After 
the  tolerance  was  entered  for  the  third  problem,  Sense  Switch  3  was 
turned  off.  When  diagonal  3  went  to  zero  and  COMPUTER  START  was 
pressed  to  continue  with  the  fourth  problem,  the  solution  began 
immediately  without  the  entry  of  a  tolerance  level,  as  Switch  3  was 
off.  Note  the  tolerance  level  in  the  warning  message  of  the  fourth 
problem  is  the  last  entry  made--that  for  problem  3»  When  this 
warning  message  appeared,  Switch  2  was  turned  off  and  calculation 
continued  to  obtain  a  solution. 


ii90£300RS 

ACCEPT  TOLERANCE  IN  FORM  *.000000  PUSH  RS 
+.OOOOOIRS 

SOLUTION 

.20000000E+Q1 

•30000000E+01 

.I4OOOOOOOE+OI 

ACCEPT  TOLERANCE  IN  FORM  +.000000  PUSH  RS 
+.000001RS 

SOLUTION 
.10000065E+01 
.19999976E+01 
•30000053E+01 
o  39999963E+01 
.I49999916E+01 
.59999990E+01 
. 70000022E+01 
. 80000012E+01 

ACCEPT  TOLERANCE  IN  FORM  +.000000  PUSH  RS 
+  .1RS 

DIAGONAL  IS  ZERO,  PUSH  START  FOR  NEXT  PROBLEM 
IS  3 

DIAGONAL  THAT  IS  LESS  THAN  OR  EQUAL  TO  TOLER 
IS  3 

TURN  SVi  2  OFF  TO  CONT.  ON  FOR  NEXT  PROBLEM 
THEN  PUSH  START 

SOLUTION 
-.67605150E+01 
.3l4l405ll5E+02 
-.70773320E+01 
.59L39liPOE+01  . 

-.69508359K+01 
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SOURCE  PROGRAM  FOR  SOLUTION  OF  SIMULTANEOUS  LINEAR  EQUATIONS 


ENTER  SOURCE  PROGRAM 
THEN  PUSH  START 

DIMENSION  X(26,27) 

30  F0RMAT(I2) 

31  FORMATTED*. 8, EUt.8,EHt.8,El]*.8,Ell*. 8) 

32  F0RMAT(//8HS0LUTI0N) 

33  FORMAT(1*1*HDIAOONAL  THAT  IS  LESS  THAN  OR  EQUAL  TO  TOLER) 
314  FORMAT(  2HIS  13) 

35  FORMATS 2HTURN  SV,  2  OFF  TO  CONT,  ON  FOR  NEXT  PROBLEM) 

36  F0RMAT(15HTH£N  PUSH  START) 

3?  FORMAT(li5HDIAOONAL  IS  ZERO,  PUSH  START  FOR  NEXT  PROBLEM) 
36  F0RMAT(l4lHACCEPT  TOLERANCE  IN  FORM  4-.  000000  PUSH  RS) 

39  FORMAT(F8.6) 
llO  F0RMAT(Ell*.8) 

10  PRINT  38 
ACCEPT  39,  TOIER 

1  READ  30,  N 
N1  -  N  4  1 

DO  2  J  •  1,  N1 
DO  2  I  -  1,  N,  5 

2  READ31,X(I,J),  X(l4l,J),  X(l42,J),  X(l43,j),  X(l4l,,J) 

DO  Ui  I  -  1,  N,  1 

DIAG  -  X(I,I) 

IF  (DIAG)  3,  20,  1) 

3  IF  (DIAG  4  TOLER)  5,  19,  19 
It  IF  (DIAG  -  TOLER)  19,  19,  5 

5  DO  6  J  -  I,  N1 

6  X(I,J)  -  X(I,J)  /  DIAG 
K  -  1 

9  IF  (K-I)  11,  13,  11 

11  FELMT  -  X(K,I) 

DO  12  J  -  I,  N1 

12  X(K, J)  -  X(K, J)  -  FELMT  *  X(I,J) 

13  K  -  K  4  1 

IF  (K-N)  9,  9,  lit 
U  CONTINUE 
J-Nl 

IF  (SENS*,  SWITCH  1)  15,  17 

15  TYPE  32 

DO  16  I  -  1,  N 
PUNCH  1*0,  X(I,J) 

16  TYPE  1*0,  X(I,J) 

GO  TO  21 

17  TYPE  32 

DO  18  1-1  ,  N 


56 


18  TYPE  1*0,  X(I,J) 

21  IF  (SENSE  SWITCH  3)  10,1 

19  TYPE  33 
TYPE  3b, I 
TYPE  35 
TYPE  36 
PAUSE 

IF  (SENSE  SWITCH  2)  21,5 

20  TYPE  37 
TYPE  3ti,I 
PAUSE 

GO  TO  2f 
END 
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